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ABSTRACT 

In this paper, we show that the pattern of microlensing hght curve anomahes 
r*^ I induced by muhiple planets are well described by the superposition of those 

of the single-planet systems where the individual planet-primary binary pairs 
act as independent lens systems. Since the outer deviation regions around 
►^ . the planetary caustics of the individual planets occur in general at different 

r~.>. ' locations, we find that the pattern of anomalies in these regions are hardly 

affected by the existence of other planet (s). This implies that even if an event 
lO ' is caused by a multiple planetary system, a simple single-planet lensing model 

C"*^ I is good enough for the description of most anomalies caused by the source's 

^— ^ . passage of the outer deviation regions. Detection of the anomalies resulting 

from the source trajectory passing both the outer deviation regions caused 
by more than two planets will provide a new channel of detecting multiple 
i-G ■ planets. 
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1 INTRODUCTION 



Since the first discovery of extra-solar planets around a pulsar PSR 1257-1-12 from the 
analysis of planetary perturbations on pulse arrival time by Wolszczan & Frail (1992) and 
c3 I Wolszczan (1994), various groups have detected planets around nearby stars using high 

precision radial velocity measurements, e.g., around 51 Pegasi (Mayor & Queloz 1995), 47 
Ursae Majoris (Butler & Marcyl996), 70 Virginis (Marcy & Butler 1996), and 16 Cygni 
B (Cochran et al. 1997). Up to now, nearly 60 extra-solar planets are known ( |http://cTa 



www. harvard. edu/planet^ ) . From close examination of the velocity residuals to the Keplerian 



fits to a subsample of 12 planet-bearing stars, Fischer et al. (2001) suggested that nearly 
1/2 of the stars may have a second companion. Currently, six systems are identified to 
have muhiple planets: around i^ And (Butler et al. 1999), HD 83443 (Mayor et al. 2001), 
HP 168443 (Marcy et al. 2001a), G7 876 (Marcy et al. 2001b), HD 82943 and HD 74156 
( |http://exoplanets.org/mult.shtm1| ). 

Although majority of the best candidate extra-solar planets were discovered by the radial 
velocity technique, they can also be detected by using other techniques. Microlensing is one of 

* clieongho@astroph. chungbuk. ac.kr 
T hyc@ns.kias. re. kr 
T jinhan@astrononiy.ohio-state.edu 
S kchang@chongju. ac.kr 



© 0000 RAS 



2 Han et al. 

the promising techniques, especially to search for planets located at Galactic-scale distances. 
Planet detections by using microlensing is possible because the event caused by a lens system 
composed of a planet (or planets) can exhibit noticeable anomalies in the lensing light curve 
when the source star passes closely to the lens caustics (Mao & Paczyhski 1991). The caustic 
refers to the source position on which the lensing-induced amplification of a point-source 
event becomes infinity. The lens system with a planet (planets) form several disconnected 
sets of caustics. Among them, one is located very close to the primary lens, "central caustic" , 
and the other(s) is (are) located away from the primary lens, "planetary caustic(s)" (Griest 
& Safizadeh 1998). The planet-induced anomalies last only a few hours to days, and thus 
it is difficult to detect them from the survey- type lensing experiments (MACHO: Alcock et 
al. 1993; OGLE: Udalski et al. 1993; EROS: Aubourg et al. 1993; MOA: Abe et al. 1997). 
However, with sufficiently frequent and accurate monitoring from follow-up observations of 
ongoing events alerted by the survey teams, one can detect the anomalies and determine 
the planet parameters of the mass ratio and the instantaneous projected separation between 
the planet and the primary lens. Currently, several groups are actually carrying out such 
observations: MPS (Rhie et al. 1999), PLANET (Albrow wt al. 1998), and MOA (Abe et 
al. 1997). 

Along with the observational efforts, planetary microlensing has also been the field of 
intensive theoretical studies. These include phenomenological studies of the lensing behav- 
iors (Wambsganss 1997; Dominik 1999; Safizadeh et al. 1999), estimation of the detection 
probability and efficiency (Gould & Loeb 1992; Bolatto & Falco 1994; Bennett & Rhie 1996; 
Peale 1997, 2001; Gaudi & Sackett 2000), development of the observational strategies for 
optimal planet detections (Griest & Safizadeh 1998; Han & Kim 2001), and various problems 
and solutions in determining planet parameters (Gaudi & Gould 1997; Gaudi 1998; Griest & 
Safizadeh 1998). However, most of these studies were focused on lens systems composed of 
only a single planet, despite the fact that our solar system is composed of multiple planets. 
Gaudi, Naber & Sackett (1998) demonstrated that if extra-solar systems are composed of 
planets having orbital separations comparable to those of Jupiter and Saturn of our solar 
system, the joint probability of two planets having projected separations within the lensing 
zone of ~ 0.6 - 1.6 of the Einstein ring radius is substantial, invoking the necessity of studies 
about the lensing properties of multiple planetary systems. 

In this paper, we investigate the properties of lensing light curve anomalies induced by 
multiple planets. Multiple planetary microlensing was first explored by Gaudi et al. (1998). 
However, since they were interested only in the anomalies occurring near the peaks of the 
light curves of high amplification events, for which the planet detection efficiencies are high 
(Griest & Safizadeh 98), they investigated anomaly patterns only in the region around 
central caustics. In this work, we extend the investigation in wider areas including the outer 
deviation regions around planetary caustics. 

The paper is organized in the following way. We begin with a brief review of the basics 
of multiple planetary lensing in § 2. In § 3, we present the anomaly patterns of an example 
lens system and the resulting light curves to illustrate the effect of multiple planets. We 
discuss about the implications of the new findings in § 4. A brief comment and conclusion 
are found in 5 5. 



2 BASICS OF MULTIPLE PLANETARY LENSING 

A lens system composed of multiple planets is described by the formalism of multiple lensing 
with very low mass-ratio companions. If a source star is lensed by A^ point-mass lenses, 
the locations of the resulting images are obtained by solving the mapping equation (lens 
equation), which is expressed in complex notations by 
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c = ^ + E-^> (1) 

where m/s represent the mass fractions of the individual lenses such that J2f itT'j = 1; ^/s 
are the positions of the lenses, C = C + ^^ ^^^ z = x + iy are the positions of the source 
and images, and z denotes the complex conjugate of z (Witt 1990). We note that all these 
lengths are normalized by the angular Einstein ring radius, which is related to the total 
mass, M, and the geometry of the lens system by 



UGM f 1 1 \i/2 

where Doi and Dos represent the distances to the lens and source from the observer, respec- 
tively. The angular separations between the microlensed images are of the order of several 
milli-arcsecs, which are too small to be resolved. However, the lensing event can be identified 
by the change of the source star flux. The amplifications of the individual images are given 
by the inverse of the determinant of the Jacobian of the lens equation evaluated at the image 
position, i.e. 

Then the total amplification is obtained by the sum of those of the individual images, i.e. 
A = Y,i ' ^i^ where Nj represents the number of images. The caustic corresponds to the 
source position on which the determinant becomes zero, i.e. |det J\ = 0. The lens equation 
describes a mapping from the lens plane onto the source plane. Then, to find image positions 
{x,y) for a given source position {(,,1]), the equation should be inverted. 

For a single lens system (A^ = 1), the lens equation can be inverted and thus is alge- 
braically solvable, yielding two solutions (and thus the same number of images). The total 
amplification is expressed in a simple form of 

A^4±^, (4) 

where u is the dimensionless lens-source separation vector normalized by 6-^. The separation 
vector is related to the single lensing parameters by 

u = (^) X + Py, (5) 

where ^e represents the time required for the source to transit 9-e (Einstein time scale), /3 is 
the closest lens-source separation in units of 9-^ (impact parameter), and to is the time at 
that moment. The unit vectors x and y are parallel with and normal to the direction of the 
relative lens-source transverse motion. 

For a multiple lens system {N > 2), on the other hand, the lens equation is non-linear 
and thus cannot be inverted. However, it can be expressed as a polynomial in z, and the 
image positions are obtained by numerically solving the polynomial. For a A^ point-mass 
lens system, the lens equation is equivalent to a (A^^ -|- l)-order polynomial in z and there 
are a maximum A^^ -|- 1 and a minimum A^ -|- 1 images and the number of images changes 
by a multiple of two as the source crosses a caustic. Due to the non-linear nature of the lens 
mapping for a multiple lens system, the lensing behavior of events caused by a multiple lens 
is not given by a simple superposition of the individual single lens events. 

Unless the planet has a separation very close to 6'e, a single-planet lens system forms a 
single central caustic and one or two planetary caustics depending on the projected sepa- 
ration between the primary and the planet. Accordingly, there exist two different types of 
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planet-induced anomalies; one affected by the planetary caustic(s) ("type I anomaly") and 
the other affected by the central caustic ("type II anomaly") (Han & Kim 2001). Compared 
to the frequency of type I anomalies, type II anomalies occur with a relatively low frequency 
due to the smaller size of the central caustic compared to that of the planetary caustic(s). 
However, the efficiency of detecting type II anomalies can be high because intensive mon- 
itoring is possible due to the known type of target events (i.e., high amplification events) 
and the predictable time of anomalies (i.e., near the peak of light curves). 

As the number of lenses increases, solving the lens equation becomes very cumbersome. 
One commonly practiced approach to obtain amplification patterns of these lens systems is 
the inverse ray-shooting technique (Schneider & Weiss 1986; Kayser, Refsdal & Stabell 1986; 
Wambsganss 1997). In this method, a large number of light rays are shoot backwards from 
the observer plane through the lens plane, and then collected (binned) in the source plane. 
Then the amplification pattern is obtained by the ratio of the surface brightness (i.e., the 
number of rays per unit area) on the source plane to that on the observer plane and the light 
curve from a particular source trajectory corresponds to the one-dimensional cut through 
the constructed amplification pattern. This method has an advantage of enabling one to 
obtain amplification patterns for an arbitrary number of lenses, but has a disadvantage of 
requiring large computation time to obtain smooth amplification patterns. 

In our analysis of multiple-planet lensing properties, we test lens systems composed 
of two Jovian- mass planets (and thus A^ = 3). This is because if extra-solar systems have 
planets with orbital separations and mass ratios similar to those of our system, other planets 
except the two Jovian planets (Jupiter and Saturn) will have mass ratios too small to affect 
the amplification pattern. In addition, the probability of more than three planets being 
simultaneously located in the lensing zone will be very small. Besides, this allows us to 
directly solve the lens equation rather than using the more time-consuming ray-shooting 
method. 



3 PROPERTIES OF ANOMALIES 

For the investigation of the properties of lensing light curve anomalies induced by multiple 
planets, we construct maps of amplification excesses. The amplification excess is defined by 

etri = , (6) 

where At^ and Aq represent the amplifications with and without planets. With the map, one 
can obtain an overview of the anomaly patterns not testing all light curves resulting from a 
large number of source trajectories. 

In the upper panel of Figure 1 , we present the constructed excess map of an example lens 
system. The lens system has two planets with mass ratios of qi = 0.003 and q2 = 0.001, which 
respectively correspond to those of Jupiter- and Saturn-mass planets around a ~ 0.3 Mq 
star. The individual planets have projected separations (normalized by ^e) of bi = 1.3 and 
62 = 0.85 and the orientation angle between the position vectors from the primary lens to 
the individual planets is = 150°. In the map, the levels of excesses are represented by 
contours, which are drawn at the levels of e = —10%, —5%, 5%, and 10% and the regions 
of positive excess are distinguished by grey scales. To show the detailed excess patterns 
in the region around the central caustic, the map of the region is expanded and presented 
separately in Figure 2. 

From the investigation of the map, one may notice an interesting point that the locations 
of the outer caustics are very similar to those of the planetary caustics of the individual 
planet-primary binary pairs. To show this similarity more closely, in the lower panel of Fig. 1 
we present the excess map constructed by superposing the excesses of the two single-planet 
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lens systems where the individual planet-primary pairs act as independent lens systems 
(binary superposition), i.e. 

EAh,i — Aq . 
A ' (7) 

=1 ^0 



1= 



where Ab,i's represent the binary- lensing amplifications of the individual planet-primary 
pairs. In Figure 3, we also present the blowup of the central region of the map. From the 
comparison of the excess patterns in the outer deviation regions, one finds that not only 
the locations of the outer caustics but also the excess patterns around them are very similar 
each other. This can be seen also from the comparison of the light curves of events caused 
by the triple lens system (solid curves in the left panels of Figure 5) to those obtained by 
the binary superposition (dotted curves in the same panels of the same figure), which are 
obtained by 



2 



Ab, = Ao{l + ebs) = V y4b,i - Ao. (8) 



=1 



From the comparison of the maps in Fig. 2 and 3, one finds that the excess patterns in the 
region around the central caustic are also similar each other, although the similarities are 
not as exact as those in the outer deviation region.^] 

To more systematically investigate how accurately the multiple-planet lensing properties 
can be represented as binary superposition, we compute the fractional deviation of the 
amplification obtained by the binary superposition from that of the exact triple lens system 
by 

^^Anp4bs_ (9) 

In Figure 4, we present the contour map of 6 for the same lens system whose excess map is 
presented in Fig. 1. We note that in the map the contours are drawn at the levels oi 6 = —1%, 
—0.5%, 0.5%, and 1%, which are much smaller than those of the contours in the map of e. In 
the right panels of Fig. 5, we also present the deviation curves of 6 for the events resulting 
from the source trajectories marked in Fig. 1. From the map and deviation curves, one finds 
that the differences occur only in small localized regions very close to the caustics. The spikes 
in the deviation curves are caused by the shift of the caustics. However, we note that the shift 
is very small (~ 10~^6'e), corresponding to the time shift of approximately half an hour for a 
typical Galactic bulge event with ^e ~ 20 - 30 days. Considering the photometric precision 
(^ 1%) and monitoring frequency (~ 3 times/night) of the current microlensing follow-up 
observations, it will be difficult to notice the caustic shifts. In addition, the modified light 
curve can be fit by a light curve with very slightly changed lensing parameters. Therefore, 
the binary superposition will be a good approximation for the description of the example 
multiple-planet system. 

Then, to what companion mass ratios is the approximation of binary superposition valid 
for the description of the lensing behavior of multiple-companion systems? To answer this 
question, we compute the fraction of the area within the Einstein ring where the binary 
superposition approximation deviates more than a threshold value ^th for lens systems with 
various combinations of mass ratios, /dev(5th)- In Figure 6, we present the contours of con- 
stant /dev in the parameter space of log qi and log g2- We note that the presented /dev are the 
values averaged for lens systems with planets having relative orientation angles in the range 
< < 27r and located in the lensing zone, i.e. 0.6 ^h <^ 1.6. The upper and lower panels 

™ For example, compared to the two disconnected sets of closed shape for the caustics constructed by the binary superposition, 
the triple lens system forms a single set of nested caustics. 
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represent the maps constructed with the threshold deviations of ^th = 0.5% and 1%, respec- 
tively. From the maps, one finds that most lens systems with heavier companion mass ratios 
of gi <^ 10~^ have /dev < 1%, implying that the lensing behaviors of multiple-planet systems 
are well approximated by the binary superposition. We note, however, that if gi ^ 0.05, 
the fraction becomes considerable [fdewi^th = 0.5%) ^ 8% and fdev{<^th = 1%) ^ 4%], im- 
plying that the binary superposition approximation will not be appropriate to describe the 

lensing behavior of systems composed of very heavy planets (with masses ^ 10 Mj)[||] and 
brown-dwarf-mass companions. 



4 IMPLICATIONS 

In the previous section, we show that the patterns of amplification excesses of a lens system 
composed of multiple planets can be well approximated by the superposition of excesses 
induced by the planets of the individual planet-primary binary pairs. This finding has the 
following implications in describing the lensing behaviors of events caused by multiple plan- 
etary systems and determining their planet parameters. 

The first implication of our finding is that even if an event is caused by a multiple-planet 
lens system, a simple single-planet lensing model will be good enough for the description of 
most type I anomalies. This is because the excess pattern in the outer deviation region caused 
by one of the planets is hardly affected by the other planet component (s). Complication 
might occur if the outer deviation regions induced by the individual planets occur at a 
similar place. However, this case is rare since the outer excess region induced by a planet- 
mass companion occupies a small area. 

However, the situation becomes complicated for the type II anomalies. This is because 
the central deviation regions of the individual planets always occur in the same region 
regardless of their orientations. As a result, the pattern of anomalies induced by one planet 
in this region can be significantly affected by the existence of other planet (s), as pointed out 
by Gaudi et al. (1998). Since one should apply complex multiple-planet lensing models for 
the description of anomalies occurred in this region, accurate determinations of the planet 
parameters from the observed light curves will be difficult, as pointed out also by Gaudi et 
al. (1998). 

5 CONCLUSION 

We have investigated the properties of anomalies in lensing light curves induced by multiple 
planets. From the analysis of the maps of excess amplification and resulting light curves, 
we find that the excess patterns of a multiple-planet lens system are well described by 
the superposition of the excesses of the single-planet systems where the individual planet- 
primary pairs act as independent lens systems. This finding has an important implication 
that a simple single-planet lensing model is good enough for the description of most type I 
anomalies even if an event is caused by a lens system composed of multiple planets. 

Another interesting point to be mentioned is that detection of the anomalies resulting 
from the source trajectory passing both the outer regions of deviation induced by more than 
two planets (e.g., the trajectory designated by a number '3' in Fig. 1) will provide a new 
channel of detecting multiple planets. The previously suggested method of detecting multiple 
planets by analyzing anomalies occurred by the source's passage of the central caustic region 
(such as the trajectory '4' in Fig. 1) has significant limitations in identifying the existence 

II For HD 168443, which was identified to have multiple planets by the radial velocity technique, one planet member has a 
mass of > 17 Mj 
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of multiple planets and determining their parameters due to the degeneracy of the resulting 
light curves from those of single-planet lensing and the complexity of multiple-planet lensing 
models. On the other hand, if a multiple-planet lens system is detected through the new 
channel, one can better determine the planet parameters because the individual anomalies 
are well approximated by relatively simple single-planet lensing models. We defer a more 
systematic analysis of this method and estimation of multiple planet detection probabilities 
for future work. 

We thank to B. S. Gaudi for making helpful comments on multiple planetary lensing. 
This work was supported by a grant (2000-015-DP0449) from the Korea Research Foundation 
(KRF). 
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Figure 1. Upper panel: The contour map of amplification excesses of a lens system composed of multiple planets. The lens 
system has two planets with the ratios of mass to that of the primary of qi = 0.003 and q2 = 0.001 and the projected separations 
from the primary of bi = 1.3 and b2 = 0.85, respectively. The locations are set so that the primary is at the center. The heavier 
planet is located on the ^ axis on the left side of the primary and the orientation angle between the position vectors to the 
individual planets (with respect to the primary) is 150°. The contours are drawn at the levels of e = —10%, —5%, 5% and 
10% and the regions of positive excesses are distinguished by grey scales. Lower panel: The excess map constructed by the 
superposition of excesses of the single-planet lens systems where the individual planet-primary pairs acts as independent lens 
systems. For both maps, the figures drawn by thick solid line represent the caustics and the straight lines with arrows are the 
source trajectories of events whose resulting light and excess curves are presented in Fig. 5. 
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Figure 2. The blowup of the central caustic region of the aniplifieation excess map presented in the upper panel of Fig. 1. 
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Figure 3. The blowup of the central caustic region of the amplification excess map presented in the lower panel of Fig. 1. 
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Figure 4. The contour map of the fractional deviation of the amplification obtained by the binary superposition approximation 
from that of the exact triple lens system. The map is for the same lens system whose excess map is presented in Fig. 1. The 
contours are drawn at the levels of 5 = —1%, —0.5%, 0.5%, and 1% and the regions of positive <5 is distinguished by grey scales. 
The lower panel is the blowup of the map in the region around the central caustic. 
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Figure 5. Left panels: The light curves of the events resulting from the source trajectories marked in Fig. 1. The number 
in each panel corresponds to the trajectory number. Right panels: The curves of the fractional deviations of the amplification 
obtained by the binary superposition approximation from that of the exact triple lens system. 
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Figure 6. The contours of constant /dovi which represents the fraction of the area within the Einstein ring where the binary 
superposition approximation deviates from the exact lensing ampUfication more than a threshold value S^^^ for the lens systems 
composed of two planets with mass ratios gi and 52 ■ The presented /dcv a^re the values averaged for lens systems with planets 
having relative orientation angle in the range < 4> < 2tt anf located in the lensing zone, i.e. 0.6 < bi < 1.6. The upper and 
lower panels are the maps constructed with the threshold deviations of (5th = 0.5% and 1%, respectively. 
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